Hyperon-nucleon single-particle potentials with low- momentum interactions 
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Single-particle potentials in Hartree-Fock approximation for different hyperon-nucleon (YN) chan- 
nels are calculated in the framework of the effective low-momentum YN interaction Vj w k- In con- 
trast to the nucleon-nucleon interaction, the available experimental data for the YN interaction are 
scarce. As a consequence no unique YN low-momentum potential Vi ow k can be predicted from the 
various bare potentials. The resulting momentum- and density-dependent single-particle potentials 
for several different bare OBE models and for chiral effective field theory are compared to each 
other. 

PACS numbers: 21.65.-f, 13.75.Ev, 21.30.Fe 



I. INTRODUCTION 

The understanding of the hyperon-nucleon (YN) in- 
teraction is essential to the physics of nuclear systems 
with strangeness and to the octet of the lightest baryons. 
Furthermore, since e.g. the core of neutron stars may con- 
tain a high fraction of hyperons a deeper knowledge of the 
YN interaction is also important for astrophysical issues. 
In recent years, an increased interest in exploring nu- 
clear systems with strangeness, especially multi-strange 
nuclear systems, becomes apparent. The YN interaction 
becomes relevant not only when the properties of simple 
hypernuclei, double strange hypernuclei but also when 
the production of hyperfragments in relativistic heavy- 
ion collisions are studied. It also strongly influences the 
composition and behavior of dense nuclear matter. For 
an recent overview of hypernuclear physics see p[. 

Unfortunately, the details of the YN interaction are 
known very poorly. On the one hand, there exist only 
a very limited amount of scattering data from which 
one could construct high-quality YN potentials based 
on meson-exchange. The existing data do not constrain 
the potentials sufficiently. On the other hand, theoret- 
ical analyses, especially for many-body systems, do not 
seem to produce unambiguous results 0. For exam- 
ple, this uncertainty is immediately demonstrated in that 
six different parameterizations of the Nijmegen YN po- 
tentials fit equally well the scattering data but produce 
very different scattering lengths, see e.g. Q. In order 
to improve the reliability of available hyperon-nucleon 
and even hyperon-hyperon potentials forthcoming exper- 
iments at the planned J-PARC and FAIR facilities are 
indispensable. Recently, first lattice QCD simulations on 
the YN interaction have been performed In [(| the 
findings on the lattice concerning some aspects of AiV 
scattering are confronted with results of the chiral effec- 
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tive field theory (%EFT ). 

For the nucleon-nucleon interaction a unique low- 
momentum effective potential is obtained 0|. This 
uniqueness is lost in the YN case due to the less con- 
straint bare potentials. Given that there are significant 
differences between various YN potentials it will be in- 
teresting to see how those differences reflect themselves 
in the hyperon single-particle potentials. 

In this work we therefore wish to compare various YN 
potentials and determine their similarities and differences 
in a dense baryonic medium. For this purpose we calcu- 
late the YN single-particle potentials and compare differ- 
ent low-energy scattering results. The potentials are ob- 
tained in a Hartree-Fock approximation with an effective 
low- momentum YN potential V\ ow k- The Kn W k is con- 
structed from several bare YN potentials [l,^. Since we 
work with isospin-symmetric nuclear matter, which con- 
tains no hyperons we can neglect the hyperon-hyperon 
(YY) interaction. This is advantageous since the YY in- 
teraction is even less known and constrained than the YN 
interaction. It has been recently shown that in the con- 
struction of the YY interaction, at leading order xEFT , 
only one additional operator, which is not present in the 
YN interaction, appears. The strength of this operator 
can be roughly estimated from existing data on double 
hypernuclei fiol ]. 

The outline of the paper is as follows: In the next Sec- 
tion we first introduce the definition and formalism of 
the effective low-momentum potential V] ow k for the YN 
interaction. The V\ ow k is obtained as a solution of an 
renormalization group equation, which needs bare YN 
potentials as initial condition. We introduce and dis- 
cuss in the following some properties of the bare YN po- 
tentials, used in this work and present some low-energy 
scattering quantities. The following section, Sec. 1111) is 
devoted to the calculation of the single-particle potentials 
in the Hartree-Fock approximation. Results for symmet- 
ric nuclear matter, the partial-wave contribution to the 
single-particle potentials and a comparison to other ap- 
proaches are shown. Finally, we end in Sec. |IV] with a 
summary and conclusion. 
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II. LOW-MOMENTUM INTERACTIONS 

We consider hyperons (Y = A, with 
strangeness S = — 1 in an infinite system of nucleons 
N composed of equal numbers of protons p and neu- 
trons n. We consider densities around nuclear satura- 
tion density po f» 0.16 fm -3 corresponding to a Fermi 
momentum kp,o = 1-35 fm -1 . The electromagnetic in- 
teraction is switched off. There exists only a very lim- 
ited amount of scattering data with which one can con- 
struct phenomenological potentials. These potentials are 
used as initial condition for solving RG equations, which 
lead to effective low- momentum YN interactions Vj ow k's. 
From these different effective Vj ow k's we finally calcu- 
late the single-particle potential for the A-hyperon. We 
start with a brief review of the construction of the low- 
momentum YN interaction. Details for several bare YN 
potentials within the RG framework used here can be 
found in Refs. @,0]. 



A. Construction of Vi ow k 

The starting point for the construction of the Vi ow k 
is the half-on-shell T-matrix, T(q', q; E y ), which is deter- 
mined by the nonrelativistic Lippmann-Schwinger equa- 
tion in momentum space. The on-shell energy is de- 
noted by E v while q' , q are relative momenta between 
a hyperon and nucleon. An effective low-momentum 
T] OW k-matrix is then obtained by introducing a cutoff 
A in the Lippmann-Schwinger kernel thus integrating 
the intermediate-state momenta up to this cutoff. At 
the same time, the bare potential in the coupled-channel 
partial-wave Lippmann-Schwinger equation are replaced 
with the corresponding low- momentum potential V[ ow k, 



£ p I ^ j2^jow fe.J/'zC^'i O^low k,zy(^ Qi ^v) 
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E y {q)-E z {l) 



(1) 



The labels y, z indicate the particle channels, e.g. y = 
YN and a, (3 denote the partial waves, e.g. a = LSJ 
where L is the angular, J the total momentum and S 
the spin. In Eq. (HJ the energies are given by 



E y (q) 



Q 

2fly ' 



(2) 



with the reduced mass \i y = MyM^ /M y and total mass 
My = My + Mjv of the hyperon My and the nucleon 
M N . 

Finally, the effective low-momentum Vl ow k is defined 
by the requirement that the T-matrices are equivalent for 
all momenta below this cutoff 



The V|o W k thus obtained is non-hermitian but never- 
theless phase-shift equivalent hermitian low-momentum 
YN interactions can be obtained. Since the low- 
momentum T-matrix, T] OW k must be cutoff- independent, 
i.e. dT]ow k/dA = an RG flow equation for Vj ow k 



dV lowk (k',k) 2V ]owk (k',A)T(A,k;A 2 



dA 



l-k 2 /A 2 



(3) 



T aa (q',q;E) 



T^ k (q',q;E) , q',q<A 



can immediately be derived. Instead of solving this flow 
equation with standard numerical methods (e.g. Runge- 
Kutta method) directly, the so-called ALS iteration 
method, pioneered by Andreozzi, Lee and Suzuki, is 
used fill \12^ . This iteration method is based on a 
similarity transformation and its solution corresponds to 
solving the flow equation. Details about the convergence 
of the ALS iteration method, applied for the coupled 
channel YN interaction, can be found in @. For the 
hyperon-nucleon interaction with strangeness S = — 1 
two different basis systems, the isospin and the particle 
basis of the bare potentials are available. We will use the 
latter. 

Furthermore, in our investigations only the diagonal 
elements of the matrices are needed. Therefore, we will 
shorten our notation further and introduce the abbrevi- 
ation Vy°t y a {q',q) — ► Vy{q) for all diagonal quantities. 



B. Bare potentials 

In order to solve the flow equation ((3]) a bare poten- 
tial as initial condition for the flow has to be chosen. In 
this work several initial YN potentials, the original Ni- 
jmegen soft core model NSC89 [3], the series of models 
NSC97a-f 3 also by the Nijmegen grou p and a recent 
model proposed by the Jiilich group [15], labeled as J04 
in the following, are used. All these models are formu- 
lated in the conventional meson-exchange (OBE) frame- 
work. They involve a set of parameters, which have to 
be determined from the available scattering data. These 
are the coupling constants of the corresponding baryon- 
baryon-meson vertices and cutoff parameters for the ver- 
tex form factors. Due to the scarce YN scattering data 
base these parameters cannot be precisely fixed in con- 
trast to the NN interaction, where a lot of scattering 
data are available. In order to construct consistently con- 
ventional OBE models for the YN interaction one usually 
assumes flavor SU(3) constraints or G-parity arguments 
on the coupling constants, and in some cases even the 
SU(6) symmetry of the quark model and adjusts their 
size by fits to NN data. The major conceptual differ- 
ence between the various conventional OBE models con- 
sists in the treatment of the scalar-meson sector, which 
plays an important role in any baryon-baryon interaction 
at intermediate ranges. In contrast to the pseudoscalar 
and vector meson sectors it is still an open issue who 
are the actual members of the lowest-lying scalar-meson 
SU(3) multiplet, what are the masses of the exchange 
particles and how, if at all, the relations for the coupling 
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constants, obtained by SU(S) flavor symmetry, should 
be applied. For example, in the older versions of the YN 
models by the Jiilich group [3, E3] a fictitious a meson 
with a mass of roughly 550 MeV, arising from correlated 
7T7T exchange was introduced. The coupling strength of 
this meson to the baryons was treated as a free parame- 
ter and finally fitted to the data. However, in the novel 
Jiilich YN potential a microscopic model of the cor- 
related 7T7T and K K exchange is established in order to fix 
the contributions in the scalar a- and vector p-channel. 
This new model incorporates also the common one-boson 
exchange parts of the lowest pseudoscalar and vector me- 
son multiplets. The corresponding coupling constants are 
determined by SU (3) flavor symmetry and the so-called 
Fj (F+D) ratios are fixed by the pseudoscalar and vector 
meson multiplets by invoking SU (6) symmetry. 

In the Nijmegen YN models, NSC89 [11], NSC97 
3 and in the recently extended soft-core model for 
strangeness S = — 2 ESC04 [3, [HI] the interaction is gen- 
erated by a genuine scalar SU (3) nonet meson exchange. 
Besides the scalar- meson nonet two additional nonets, 
the pseudoscalar and vector SU(3) flavor nonets are con- 
sidered in all Nijmegen models. In addition, Pomeron ex- 
change is also included, which provides additional short- 
range repulsion. Nevertheless, there are a few concep- 
tual differences in the various Nijmegen models. In the 
NSC97 models the strength parameter for the spin-spin 
interaction the magnetic Fj (F + D) ratio is left as an 
open parameter and takes six different values in between 
a range of 0.4447 to 0.3647 for the six different mod- 
els NSC97a-f. In the original Nijmegen SC89 model this 
parameter is constrained by weak-decay data. Further- 
more, the NSC97 models include additional SU(3) flavor 
breaking, which is based on the so-called 3 Po model [2C| . 

The predictions of the above mentioned models can 
be compared with an another approach, the so-called 
chiral effective field theory (%EFT ) to nuclear interac- 
tions, which is based on chiral perturbation theory (for 
recent reviews see e.g. (2ll.[2^.[23T|). The major benefit of 
the xEFT is the underlying power counting scheme, pro- 
posed by Weinberg [2J, [23], that allows to improve the 
calculations systematically by going to higher orders in 
the expansion. In addition, higher two- and three-body 
forces can be derived consistently in this framework. Fur- 
thermore, the effective potential is explicitly energy in- 
dependent in contrast to the original Weinberg scheme. 

Within xEFT the NN interaction has been analyzed 
recently to a high precision (N 3 LO) jH]. To leading 
order (LO) the NN potential is composed of pion ex- 
changes and a series of contact interactions with an in- 
creasing number of derivatives which parameterize the 
singular short-range part of the NN force. In order to 
remove the high-energy components of the baryonic and 
pseudoscalar meson fields a cutoff A dependent regulator 
function in the Lippmann-Schwinger (LS) equation is in- 
troduced. With this regularized LS equation observable 
quantities can be calculated. The cutoff range is limited 
from below by the mass of the pseudoscalar exchange 



mesons. Note, that in conventional meson-exchange 
models the LS equation is not regularized and conver- 
gence is achieved by introducing form factors with cor- 
responding cutoff masses for each meson-baryon-baryon 
vertices. 

So far, the YN interaction has not been investigated 
in the context of the xEFT as extensively as the NN 
interaction. A recent application to the YN interaction 
by the Jiilich group can be found e.g. in [22j]. Analo- 
gous to the NN case, the YN potential, obtained in LO 
XEFT , and consists of four-baryon contact terms and 
pseudoscalar meson (Goldstone boson) exchanges, which 
are all related by SU (3) / symmetry. For the YN interac- 
tion typical values for the cutoff lie in the range between 
550 and 700 MeV (see e.g. (H). At LO xEFT and for 
a fixed cutoff A and pseudoscalar Fj (F + D) ratio there 
are five free parameters. The remaining interaction in the 
other YN channels are then determined by SU{2>) / sym- 
metry. A next-to-leading order (NLO) xEFT analysis 
of the YN scattering and of the hyperon mass shifts in 
nuclear matter was performed in [28|. However, in this 
analysis the pseudoscalar meson-exchange contributions 
were not taken into account explicitly but the YN scat- 
tering data could be described successfully for laboratory 
momenta below 200 MeV using 12 free parameters. One 
ambiguity in this approach is the value of the ry coupling 
which is identified with the octet 773 meson coupling and 
not with the physical 77 meson. The influence of this am- 
biguity on the data description can be disregarded [29| . 

Since there are scarce YN scattering data available, 
it has not been possible yet to determine uniquely the 
spin structure of the YN interaction. Nevertheless, all 
of the above mentioned OBE models are consistent with 
the measured YN scattering observables. In addition, 
all of these potentials include the AA^ — SAT conversion 
process. 



C. Low-energy scattering 

In order to obtain further insight into the separation 
of scales for the evolution of the low-momentum V] OW k 
we investigate its cutoff dependence. A common fea- 
ture of all YN potentials is the long-range one-pion 
exchange (OPE) tail. In general, the RG decimation 
eliminates the short-distance part of the bare potential 
and preserves the model-independent impact of the high- 
momentum components on low-momentum observables. 
In this sense, the ambiguities associated with the unre- 
solved short-distance parts of the interaction disappear 
and a universal low-momentum YN interaction V\ ow k 
can be constructed from phase shift equivalent bare YN 
potentials. 

The mentioned hierarchy of scales can be seen e.g. in 
the E~n channel, see Fig. [TJ The Vj ow k matrix elements 
for vanishing momenta are shown as a function of the cut- 
off A for the 1 S'o partial wave. When A is decreased, the 
resulting VJ ow k becomes more and more attractive. For 
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1 5'o and a cutoff A ~ 500 — 250 MeV V\ ow k becomes cut- 
off independent. Decreasing the cutoff further below the 
2ir exchange threshold, which corresponds to a k « 280 
MeV, the cutoff insensitivity disappears since the pion 
contributions are finally integrated out. 
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FIG. 1: Vio W k(0) i n 1, S'o partial wave for various bare po- 
tentials as a function of the cutoff A in the E~n channel. 
Prediction from effective range theory (lines) are added. 

In the opposite direction, i.e. for A — > oo no fluctua- 
tions have been integrated and V] OW k tends to the bare 
potential. 

The limit A — > of Mow k should yield the scattering 
length. In the limit of small cutoffs an analytic solution 
obtained in the framework of the effective theory, see [3p| . 
is given by the expression 



Vy(0) 



2^-2± 
a 7T 



for 



A^0 



(4) 



where we have simplified our notation in an obvious man- 
ner. Here, the scattering length ag is needed as input, 
which we have calculated in the standard effective range 
approximation directly from the T-matrix for the So 
channel from the Vj ow k- In this approximation the T- 
matrix for q < A can be expanded as 



q cot do = — 



1 



2fi y T y (q, q; q 2 ) 



1 

a 



1 



-joq 



(5) 



where ro is the effective range. The results for all differ- 
ent YN flavor channels and for all bare OBE potentials 
and the xEFT potential used in this work with cutoffs 
between 550 and 700 MeV are listed in Tab. UJ for the 
scattering length ao in units of fm and in Tab. ITU for the 
effective range ro also in units of fm. 

As can be seen from Fig. Q] there is good agreement 
for small cutoff values A between the analytical expan- 
sion and the full V[ ow k solution obtained from the flow 
equation. 





Ap 


An 


E°p 


E°n 


E+p 


E+n 


E"p 


E"n 


NSC97a 


-0.71 


-0.76 


-2.46 


-1.74 


-6.06 


-0.04 


0.41 


-6.13 


NSC97b 


-0.90 


-0.96 


-2.47 


-1.72 


-5.98 


-0.04 


0.41 


-6.06 


NSC97c 


-1.20 


-1.28 


-2.41 


-1.70 


-5.90 


-0.03 


0.41 


-5.98 


NSC97d 


-1.70 


-1.82 


-2.38 


-1.68 


-5.82 


-0.03 


0.41 


-5.89 


NSC97e 


-2.10 


-2.24 


-2.38 


-1.68 


-5.82 


-0.03 


0.41 


-5.90 


NSC97f 


-2.51 


-2.68 


-2.45 


-1.74 


-6.07 


-0.05 


0.42 


-6.16 


NSC89 


-2.70 


-2.72 


-2.12 


-1.57 


-4.79 


-0.09 


0.23 


-4.85 


J04 


-2.14 


-2.11 


-2.24 


-1.63 


-4.68 


-0.18 


0.04 


-4.75 


XEFT550 


-1.80 


-1.79 


-1.76 


-1.15 


-3.82 


0.12 


0.31 


-3.88 


XEFT600 


-1.80 


-1.80 


-1.25 


-0.92 


-2.70 


0.10 


0.20 


-2.72 


XEFT650 


-1.80 


-1.80 


-1.43 


-1.02 


-3.06 


0.09 


0.21 


-3.10 


XEFT700 


-1.80 


-1.80 


-1.50 


-1.07 


-3.19 


0.06 


0.20 


-3.24 



TABLE I: Scattering lengths ao of Vi ow k for different flavor 
channels in units of fm for the 1 So partial wave. 





Ap 


An 




E°n 


E+p 


E+n 


E-p 


E"n 


NSC97a 


5.87 


6.12 


4.58 


0.60 


3.28 


-6602 


24.8 


3.27 


NSC97b 


4.93 


5.10 


4.68 


0.59 


3.29 


-8491 


25.0 


3.28 


NSC97c 


4.11 


4.23 


4.79 


0.57 


3.30 


-10670 


25.4 


3.29 


NSC97d 


3.46 


3.53 


4.91 


0.54 


3.30 


-17115 


25.4 


3.29 


NSC97e 


3.19 


3.24 


4.90 


0.52 


3.29 


-17326 


25.2 


3.29 


NSC97f 


3.03 


3.09 


4.60 


0.51 


3.25 


-6341 


24.1 


3.24 


NSC89 


2.86 


2.98 


5.76 


0.74 


3.35 


-1478 


58.0 


3.33 


J04 


2.93 


3.09 


3.76 


1.04 


3.32 


-329 


1232.0 


3.30 


XEFT550 


1.73 


1.84 


6.10 


-2.96 


2.70 


-825 


34.1 


2.68 


XEFT600 


1.77 


1.88 


5.32 


-2.12 


3.40 


-780 


10.2 


3.39 


XEFT650 


1.75 


1.86 


5.10 


-2.28 


3.08 


-1210 


27.6 


3.05 


XEFT700 


1.74 


1.86 


4.91 


-2.17 


2.97 


-2450 


34.8 


2.95 



TABLE II: Effective range ro. Labeling is the same as in 
Tib. II 



Unfortunately, no general quantitative conclusion can 
be drawn from Tab. UJ and Tab. ILTIdue to the bad exper- 
imental situation for the YN data. The YN interaction 
is yet largely unknown. However, agreement of the scat- 
tering length of all NSC97 potentials except for the Ap 
and An channels is found. This deviation is related to 
the different fits of the magnetic F/(F + D) ratio in the 
Nijmegen potentials Q. The remaining two potentials, 
NSC89 and J04, have different but comparable values to 
those of the NSC97 ones. Unfortunately, the difference 
between these potentials and the xEFT potential is large. 

Fig. [2] shows the same matrix element as in Fig. [T] but 
for the 3 S'i partial wave. Unlike to the 1 So channel, Vj ow k 
for the Si channel is still cutoff dependent. This can 
be traced back to the different short-range behavior of 
the two channels. In the 1 S channel the potential has 
an strongly repulsive core while in the 3 Si channel the 
short-distance part is strongly attractive. Hence, dur- 
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ing the RG decimation towards smaller cutoff values, the 
potential gets more and more attractive. 
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FIG. 2: V^ w k (0) similar to Fig. Q] for the 3 Si channel. 



III. SINGLE-PARTICLE POTENTIALS 

Generally, the single-particle potential U (p) is defined 
as the diagonal part in spin and isospin space of the 
proper self-energy for the single-particle Green's func- 
tion. In the Hartree-Fock approximation for a uniform 
system it represents the first-order the interaction energy 
of a particle with incoming momentum p and given spin 
and isospin with the filled Fermi sea. For the YN interac- 
tion the hyperon-nucleon single-particle potential Uy(p) 
describes the behavior of the incoming hyperon Y with 
momentum p in the nuclear medium, i.e. its interaction 
with a filled Fermi sea of nucleons. The ground-state 
energy of nuclear matter in first approximation is then 
obtained by integrating U(p) over the incoming momen- 
tum, spin and isospin and adding the kinetic energy. Pic- 
torially, the ground-state energy is represented by closed 
Goldstone diagrams. By cutting one line, symbolizing 
the hyperon propagator, in each corresponding Goldstone 
diagram the Uy(p) is obtained. 

In the following, the calculation of the momentum- and 
density-dependent Uy{p) in the Hartree-Fock approxima- 
tion is presented for symmetric nuclear matter. 



A. Hartree-Fock approximation 

The single-particle potential Uy{p) for a hyperon with 
momentum p = \p\ is obtained from the diagonal el- 
ements of the low-momentum potential matrix, V™(q), 
where the labeling introduced in Eq. (fl]) has been used. 



In a plane- wave basis it splits into two terms, the (direct) 
Hartree- and the (exchange) Fock-term [3l| . 

KM = WLect + (-!) L+S ^)L chango ■ (6) 

Since the low-momentum potential is given in the partial- 
wave basis we perform a change of basis with the result 



Uy(p) = 



E 

a, Ms, Ml 



rj2 fdtP 2 LML (t) fd qq 2 v«( q ) , (?) 



where q — (Myk — M^p)jM y denotes the relative mo- 
mentum between the nucleon k and hyperon momentum 
p. The Plm l (*) are the associated Legendre polynomials 
of the first kind with the argument t = p-q, where the hat 
labels an unit vector. The quantity f a is the short-hand 
notation for the square of the Clebsch-Gordon coefficients 



r 



r 


S 






M s 


Mj j 



2L + l(L- M L )\ 
2tt (L + M l )\ 



(8) 



and the Mj, Ml and Ms are the corresponding projec- 
tions of the total J, angular momentum L and spin S, 
respectively. 

The integration boundaries of Uy(p) in Eq. ([7]) are 
derived in the following: For vanishing temperature the 
momentum of the nucleon k is restricted by the Fermi 
momentum kp, i.e. < k < kp, and kp is directly 
related to the proton or neutron density pjv via 



k F = 3tt 2 XnPb 



(9) 



Here xjv = Pn/pb describes the ratio between the proton 
or neutron and total baryon density. The inequality < 
k < kp can be reformulated as 

M 2 q 2 + M%p 2 + 2M y M N qpt - M Y kp < 

which has the solution for the relative momentum q 

q~(k F ,p,t) <q< q + (k F ,p,t) 

with the definitions 



q ± {k Fl p,t) 



N 



M 
~M V 



p-t± 



M- 
~M 



N 



(10) 

Since the relative momentum q is a real quantity this 
further constrains the integration variable t to 



t>\ 1- 



/ My kp 

fe7 



(11) 



which is only valid if the hyperon momentum is p > 
jj^-kp. In this case this finally determines the integra- 
tion limits as 



/ My kp 

q m in =q~{k F ,p,t) 



q + {kF :Pl t) (12) 



because the modulus of t is always smaller or equal one. 



For the case that the hyperon momenta p < jf-kp the 
functions q ± (k F ,p, t) are always real which then yield the 
integration limits 



Qirtin 



Qmax = q + (k F ,p,t) 



(13) 



U A (p) [MeV] 



In Fig. [3] the integration limit functions q ± are shown 
as a function of t for three different choices of the hyperon 
momentum (p < jj^-kp,p > jj^kp,p = j^kp) and a 
fixed kp. 



t q±(t) 



a 8 



NfcS o 

* + + tC 1 



FIG. 3: q ± (kp,p,t) as function of t for different choices of p 
and fixed kp. 

Thus, the integration limits are known analytically and 
finally the Uy{p) are calculated via Eq. (7} numerically 
with standard integration methods. 



B. Symmetric nuclear matter 

For symmetric nuclear matter the ratio in Eq. © has 
to be fixed to x n = 1/2. As an example, the numerical 
solution of Eq. (7]) for the full momentum and density de- 
pendent single-particle potential of the A hyperon with 
momenta up to 500 MeV and nuclear densities up to 6po 
is shown in Fig. [4] where as bare potential for the underly- 
ing V| OW k calculation the NSC97f model of the Nijmegen 
group has been used. One sees that with increasing den- 
sity, the momentum dependence becomes stronger, indi- 
cating that the effective mass becomes lower for larger 
densities. 

Similarly, Fig. [5] shows the full momentum and den- 
sity dependence of the E~ potential in symmetric nuclear 
matter based on the NSC97f potential. Here, the slope 
of the momentum dependence is less pronounced, which 




p[Mevf 



50$ 



FIG. 4: Momentum and density dependence of Ua(p) for sym- 
metric nuclear matter. As bare potential the NSC97f was 
used. 
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FIG. 5: Similar to Fig.Gflfor U s -(p), 



leads to a weaker density dependent effective mass. How- 
ever, unlike to the A case, the curvature becomes nega- 
tive at higher densities leading to an larger effective mass 
than the bare mass. 
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FIG. 7: Momentum dependence of U\ (p) at saturation den- 
sity in symmetric nuclear matter. 



The density dependence for several A potentials at rest 
(i.e. p = 0) in symmetric nuclear matter is compared in 
Fig. [6j The square represents the generally excepted em- 
pirical potential depth of U\(p = 0) ps —30 MeV. This 
value has been recently confirmed by an analysis of the 
(ir~ , K + ) inclusive spectra on various target nuclei as 
best fits in a framework of a distorted-wave approxima- 
tion j32j. While most potentials can reproduce this value, 
the Jiilich potential (J04) yields a stronger binding while 
the old Nijmegen potential (NSC89) underestimates the 
binding. 

With the exception of the J04 and NSC89 potentials, 
all others agree up to the saturation density. However, 
with increasing density, the differences between these po- 
tentials grow, leading to different bindings at rest. This 
will have consequences for the predictions of the A hy- 
peron concentration in dense nuclear matter. In particu- 
lar, this will affect the maximum mass of neutron stars. 
It is interesting to observe that even the Nijmegen po- 
tentials NSC97a-f are different from each other at higher 
densities, since the only difference between the bare po- 
tentials NSC97a-f is the F/(F + D) ratio. 

In the past, the potentials NSC89 [H and NSC97a,f 
[3| have also been used as a basis for a single-particle 
potential calculation in the G-matrix formalism. These 
G-matrix calculations yield a more attractive A potential. 
For example, at saturation density a potential depth of 
-29.8 MeV is found for the NSC89 potential, the NSC97a 
gives -39.7 MeV, and the NSC97f -36.6 MeV. On the 
other hand, a comparison with another G-matrix calcu- 
lation Q, which uses a different prescription for interme- 
diate state spectra, yields similar results to ours. 

In Fig.[7|the momentum dependence of the A potential 



at saturation density for various YN potentials is shown. 
While all potentials increase with increasing momentum, 
the slopes deviate of each other. Similar differences in 
the momentum behavior of the single^particle potentials 
are also seen in other works, cf. e.g. [H, HH- 
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3 S 1 


'Pi 


3 Po 


3 Pi 


3 P 2 


3 D 1 


Ua 


NSC97a 


-4.86 


-27.79 


1.70 


-0.10 


2.10 


-2.03 


-0.09 


-32.12 


NSC97b 


-6.69 


-27.40 


1.86 


0.05 


2.53 


-1.87 


-0.09 


-32.72 


NSC97c 


-9.06 


-27.54 


1.96 


0.36 


2.84 


-1.72 


-0.09 


-34.42 


NSC97d 


-12.14 


-26.05 


2.22 


0.64 


3.54 


-1.33 


-0.08 


-34.46 


NSC97e 


-13.92 


-24.43 


2.43 


0.75 


4.09 


-1.03 


-0.07 


-33.50 


NSC97f 


-15.37 


-20.85 


2.85 


0.68 


5.09 


-0.47 


-0.05 


-29.49 


NSC89 


-15.73 


4.52 


2.00 


0.52 


2.55 


-3.46 


-0.07 


-10.84 


J04 


-9.55 


-35.18 


-0.15 


-0.70 


0.58 


-3.17 


-1.31 


-50.28 


XEFT550 


-11.11 


-15.46 


1.50 


-1.69 


3.17 


-0.07 


-3.14 


-27.14 


XEFT600 


-12.29 


-11.39 


1.50 


-1.73 


3.17 


-0.07 


-6.14 


-27.37 


XEFT650 


-11.99 


-6.70 


1.50 


-1.77 


3.17 


-0.07 


-9.90 


-26.27 


XEFT700 


-11.91 


-1.77 


1.50 


-1.81 


3.17 


-0.08 


-13.84 


-25.35 



TABLE III: Partial wave contributions to the A potential 
Ua(p = 0) at pB = po in symmetric nuclear matter. 



In addition, Eq. (7| cannot only be used for the calcu- 
lation of Uy(p), but also to extract the individual partial- 
wave contributions to the total potential. These contri- 
butions are obtained by neglecting the summation over 
the LSJ quantum numbers in Eq. {£j|, which we label in 
the following as U Y { 2S+1 Lj). In Tab. IHD and Tab. El 
the resulting partial-wave contributions to (7a and Us-, 
respectively, are listed for several YN interactions at van- 
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ishing momenta at saturation density. 





^0 


3 5i 


'Pi 


3 P> 


3 Pi 


3 P 2 


3 Di 


U s - 


NSC97a 


3.51 


-4.87 


-2.16 


0.59 
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-2.41 


-0.01 


-4.73 


NSC97b 


3.58 
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-6.50 
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0.00 


-5.86 


NSC97d 


3.50 


-6.08 


-2.02 


0.71 


1.70 


-1.92 
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0.28 
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J04 
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-3.37 
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-15.13 


XEFT550 


2.28 
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-2.73 


14.11 
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-3.70 


66.26 


1.50 


-0.28 


0.06 


-0.01 


-5.36 


56.89 


XEFT650 


-2.72 


42.41 


1.50 


-0.35 


0.01 


-0.01 


-8.60 


30.38 


XEFT700 


-2.93 


39.93 


1.50 


-0.41 


-0.04 


-0.02 


-11.60 


24.68 



TABLE IV: Partial-wave contributions to the E single- 
particle potential U s - (p = 0) at ps = pa. 




P [fnr 3 ] 

FIG. 8: Density dependence of Ua(p = 0) for symmetric nu- 
clear matter, p = 2k%/3TY 2 . Full lines are from [33] and dashed 
lines represent xEFT for various regulator cutoffs. 



In these tables the partial waves up to L = 2 are shown 
and the last column contains the sum up to L = 5. As 
expected, the influence of the S- and D-waves is the most 
dominant one. One can see that the combination of the 
coupled 3 Si and 3 Di channels provides most of the at- 
traction in the majority of the A potentials. 

These tables also illustrate the different contributions 
to the hyperon potentials originating from the central, 
spin-spin, spin-orbit parts of the YN interaction. These 
partial waves can than be use to determine the size of 
these contributions, which could than be compared to 
other results such as 0, [35|, [H, H?}- Furthermore, one 
can recognize that a change in the F j (F + D) ratio from 
the different bare NSC97a-f potentials affects U\ stronger 
than the Us,. 

Another interesting feature is that the 
XEFT successfully reproduces the potential depth at 
saturation density. For these densities, the xEFT agrees 
well with the Nijmegen NSC97a-f potentials. Fig. 
shows a comparison of the U\(jp = 0) density depen- 
dence, obtained with the xEFT , with results from 
[37I . Perfect agreement for the U\(p = 0) is evident 
and the independence of the xEFT potential on the 
regulator cutoff is also seen. This suggests that the two 
approaches in Refs. and to construct an xEFT , 
are closely related. Furthermore, xEFT in leading order 
can already produce a reasonable A7V potential. 

Fig. [9] shows the density dependence for several E~ 
potentials at rest in symmetric nuclear matter similar to 
the previous figure. The other members of the £ triplet, 
S + and E°, exhibit an almost identical behavior. A small 
differences compared to the S~ case is found to be due 
to a small difference in their masses. Therefore, in the 
remaining section we discuss only the E~ potential. For 
U s - no density range is found where all or even most 
potentials agree. However, the difference between the 



350 r 

: NSC97a 
300 -NSC97C 

: NSC97f 



20Q XEFT600 - - - : 

> 

CD 

11 150 

o y 
A 100 : y 

u 

=> 50 : y - 

o _'_ 

-50 : 

-100 r " 

1 2 3 4 5 6 

PEstPo] 

FIG. 9: Similar to Fig.[6]for the (7 B - (p = 0). 

NSC97a-f potentials is not significant and is the same 
over the entire range of densities shown. This further 
confirms that the influence of the F j (F + D) ratio on the 
T,N interaction is less important than on A7V. Due to 
experimental uncertainties in case of the £ potential, no 
generally accepted empirical point can be used as a refer- 
ence. On the one hand, recent results based on the 
distorted-wave impulse approximation, yield a repulsive 
potential of the order of 100 MeV. On the other hand, 
the analysis of the same data by [3^| in a semiclassical 
distorted-wave model and an analysis by within a 
distorted-wave impulse approximation with a local opti- 
mal Fermi averaged T-matrix find a less repulsive poten- 
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tial. In addition, there exists a bound state of %lie [4lj . 
which definitely requires an attractive potential. Thus, 
for this case no conclusive answer on the theoretical as 
well as experimental side can be given. 

Compared with a G-matrix calculation a stronger 
binding for the S single particle potential is found. In 
particular, using the NSC89 interaction a binding of 
-15.3 MeV [13 is found while for the NSC97a poten- 
tial -29.7 MeV and for the NSC97f potential -25.5 MeV 
are reported . In order to understand the origin of such 
a significant difference the individual partial- wave contri- 
butions to the potential by [33| are compared with each 
other. The 1 5o channel contributions are approximately 
the same while those for the 3 Si channel are significantly 
different. 

This difference in the 3 Si channel is present for both 
A and E~ potentials and is the results of the difference 
in the treatment of the 3 Si AN — channel. Since the 
effective interactions are constructed from Vj ow k and the 
G-matrix formalisms, respectively, the difference comes 
from the treatment of the attractive part of the bare po- 
tentials above the cutoff. This is similar to the point 
made in Fig. [2] where a cutoff dependence is visible: by 
lowering the cutoff more and more 'attraction' is effec- 
tively added to the interaction. It is also interesting to 
note, that the effective potentials constructed in the G- 
matrix calculations, NSC89 and NSC97a,f depend on the 
underlying bare potentials in a similar way as the poten- 
tials shown here. This is another sign that the uncertain- 
ties are inherent in the underlying potentials. 
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FIG. 10: Similar to Fig. for U s - (p). 

In Fig. [10] the momentum dependence of the E/j;- at 
saturation density for various YN potentials is displayed 
and illustrates how strongly the results depend on the 
parameterization of the YN interaction and reiterates 
that the E7V interaction is very poorly constrained. 



> 

CD 



100 
90 
80 
70 
60 
50 
40 
30 
20 
10 


-10 
-20 
-30 



- 






^EFTGOO] 




XEFTQ50: 






J* _ - 


XEFT550 


- W(k,) 









0.05 0.1 0.15 0.2 0.25 

P [fm- 3 ] 



FIG. 11: Density dependence of Us (p = 0) for symmetric nu- 
clear matter. The full line comes from [35[ and the dashed 
lines represent xEFT results for various values of the regula- 
tor cutoff. 



Fig. [UJ shows the density dependence of the real part 
Us{p = 0) and the imaginary WyAp = 0) of an opti- 
cal potential calculation of Ref. [35| together with the 
results obtained from xEFT . The most interesting fea- 
ture here is that all potentials are positive and grow with 
increasing density in contrast to other potentials. How- 
ever, unlike in the case of the A potential, f/ s - depends 
on the regulator cutoff and only x^FT with a cutoff of 
600 MeV agrees with the results of Ref. [35[ quantita- 
tively. As already mentioned earlier, the repulsive E~ 
potential, which grows with density, has been suggested 
by Saha et al. [32] by means of an analysis of (tt~ , K + ) 
inclusive spectra. 

Recently, a calculation of the binding energy of 
the A hyperon in nuclear matter within a Dirac- 
Brueckner-Hartree-Fock framework was performed us- 
ing the most recent Julich meson exchange YN po- 
tential [i^] . The reported values of the A potential 
of -51.27 MeV (-47.4 MeV) in Brueckner-Hartree- 
Fock (Dirac-Brueckner-Hartree-Fock) framework agree 
well with our prediction of —50.28 MeV. 



IV. SUMMARY AND CONCLUSIONS 

In this paper an application of recently constructed 
YN V\ ow k potentials [3, 0| is presented. The potentials 
were constructed in a RG formalism and applied to nu- 
clear matter. The calculation of the single-particle po- 
tential of the A and £ hyperon was performed in the 
Hartree-Fock approximation. Since the Vj ow k is an low- 
momentum interaction, it can be directly used in the 
considered range of momenta. 

For comparison with other approaches, scattering 
lengths and the single-particle potential of several bare 
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potentials were calculated. In contrast to the NN case, 
it is not possible to obtain a precise fit of the YN po- 
tential to all partial waves, due to the incomplete YN 
scattering data base. This deficiency is clearly visible in 
the results obtained. 

A few interesting observations are in order, however. 
The V[ ow kfor the £~V 1 S , channel is cutoff indepen- 
dent in between momenta of 250 and 500 MeV, but no 
such independence was found in the 3 Si channel due to 
a strongly attractive short-range part. The results for 
the A potential are more satisfying since most of the 
potentials agree up to saturation density and reproduce 
the empirical point. An interesting observation is that 
%EFT in leading order lies among these potentials. How- 
ever, for the E~ potential no such agreement exists. Un- 
fortunately, the present experimental situation for S~ 
does not allow for more stringent conclusions. 

Finally, the results are also compared with other ap- 
proaches, [1, H3, S|- Most of the results of a G-matrix 
calculations yield a stronger binding than we find. The 
difference comes from the way in which the attraction, es- 
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